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We present theoretical predictions of spintronic transport phenomena that should be observable 
in ferromagnetic Fe nanocontacts bridged by chains of Fe atoms. We develop appropriate model 
Hamiltonians based on semi-empirical considerations and the known electronic structure of bulk Fe 
derived from ab initio density functional calculations. Our model is shown to provide a satisfactory 
description of the surface properties of Fe nano-clusters as well as bulk properties. Lippmann- 
Schwinger and Green's function techniques are used together with Landauer theory to predict the 
current, magneto-resistance, and spin polarization of the current in Fe nanocontacts bridged by 
atomic chains under applied bias. Unusual device characteristics are predicted including negative 
magneto-resistance and spin polarization of the current, as well as spin polarization of the current 
for anti-parallel magnetization of the Fe nanocontacts under moderate applied bias. We explore the 
effects that stretching the atomic chain has on the magneto-resistance and spin polarization and 
predict a cross-over regime in which the spin polarization of the current for parallel magnetization of 
the contacts switches from negative to positive. We find resonant transmission due to dangling bond 
formation on tip atoms as the chain is stretched through its breaking point to play an important 
role in spin-dependent transport in this regime. The physical mechanisms underlying the predicted 
phenomena are discussed. 

PACS numbers: 75.47.-m,73.63.-b,72.25.-b,72.10.-d 



I. INTRODUCTION 

Spin-dependent transport (SDT) phenomena occur 
when a bias voltage is applied across a junction between 
materials one or both of which are magnetic i They in- 
clude partial or complete spin-polarization of the electric 
current and changes in the electrical resistance when the 
magnetization direction of one of the magnetic compo- 
nents of the system is reversed through the application 
of a magnetic field. The latter effect is usually referred 
to as "junction magneto-resistance" or "giant magneto- 
resistance" (GMR). SDT has been observed at interfaces 
between ferromagnetic metals and superconductors^, be- 
tween ferromagnetic and normal metals^, between fer- 
romagnetic metals separated by thin insulating films 4 , 
and more recently between magnetic semiconductors and 
nonmagnetic semiconductors^^ 

At the present time spin-dependent transport through 
nanoscale junctions is attracting increasing attention. 
SDT through molecules bridging nanoscale magnetic con- 
tacts has been investigated theoretically2i2i2ii£ and rele- 
vant experiments have recently been reportediii*i2ii2ii4 
SDT through chains of atoms or single atoms con- 
necting pairs of nickel and cobalt nanocontacts has 
also been investigated theoreticallyii^4iS*iii2iiS Such sys- 
tems have been realized experimentally using break- 
junction22*2ii22i- 23 ' 24 and electrochemical techniques 
and transport measurements on them have been carried 
outiSSiSiiS^iS&S 4 * 2 ^ SDT measurements have also been car- 
ried out on Fe atomic contacts and magnetoresistances an 
order of magnitude smaller than for the Ni and Co sys- 
tems were founds This difference was explained in terms 
of the smaller ratio between the spin up and spin down 
densities of states at the Fermi level in bulk Fe than in Ni 



or Coi2& Very recently, a theoretical study elucidating the 
systematics of SDT through a variety of nanocontacts, in- 
cluding Fe nanocontacts connected by atomic chains, has 
been reported^ The calculations were based on a semi- 
empirical tight-binding model and were carried out in the 
limiting case of infinitesimal applied bias voltage. The 
structures studied were periodic arrays of nanoconstric- 
tions with all of the Fe atoms occupying sites of a bulk 
crystal lattice and enough vacant sites being included in 
the structure to isolate the Fe chains from each other. 

In this article we explore SDT through Fe atomic 
chains connecting Fe nanocontacts theoretically in some 
important regimes that were not considered in the previ- 
ous work: We develop a model applicable to more general 
atomic geometries; thus we are able to examine the ef- 
fects on SDT of stretching the atomic chain that bridges 
the Fe nanocontacts as occurs in break-junction experi- 
ments. We predict pronounced resonant spin-dependent 
transport phenomena due to dangling bonds that form 
on the tip atoms as the atomic chain parts. We also con- 
sider the application of finite bias across the junction of 
the Fe nanocontacts and predict other unusual SDT phe- 
nomena including negative junction magneto-resistance 
and negative spin-polarization of the current. We define 
these quantities as follows: 

The junction magneto-resistance (JMR) is defined as 

JMR = [ Ipar - Iantl \ (1) 

where I par (I an u) is the electric current flowing between 
the Fe nanocontacts when their magnetizations are par- 
allel (antiparallel). 
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The spin-polarization of the current (SP) is defined as 

gp _ {lup ~ I down) ^ 
2\^up "T" Idown) 

Note that we refer to the majority spins in eac/i contact as 
"spin up" and the minority spins as "spin down" . Thus 
if an electron is transmitted between contacts with anti- 
parallel magnetizations without changing its spin orien- 
tation, the transition will be referred to as a spin up — * 
spin down or spin down — > spin up transition. In equa- 
tion J2J I up {I down) is the spin up (spin down) current 
with up and down defined as for the nanocontact that 
is the electron drain electrode for the system. The total 
current is defined as / = I up + Idown for both parallel and 
antiparallel magnetizations. 

We find that, interestingly, in equation J3J) the cur- 
rent for the spin with the larger bulk density of states is 
not always larger than the current for the spin with the 
smaller density of states. Thus the sign of the spin polar- 
ization of the current cannot be predicted from a knowl- 
edge of the spin up and down densities of states alone, 
even when the magnetizations of the two nanocontacts 
are parallel. We also show that the spin polarization of 
the current need not vanish even for the case of anti- 
parallel magnetizations of the two nanocontacts because 
of a symmetry breaking that occurs in the system un- 
der finite applied bias. Furthermore we predict that the 
sign of the the junction magnetoresistance should change 
when the atomic chain connecting the nanocontacts is 
stretched and also if the applied bias voltage is increased 
sufficiently. 

Our SDT calculations are based on Landauer theorjiSi 
and Lippmann-Schwinger and Green's function tech- 
niques. Transport calculations require a knowledge of 
the underlying electronic structure. Semi-empirical tight- 
binding models of electronic structure have been used 
successfully in modeling SDT in thin film structures in- 
volving Fe and other magnetic metals together with insu- 
lating or vacuum tunnel barriers or non-magnetic metal 
spacers j 28 i 29 ' 30 i 31 i 32 and also SDT in magnetic metal 
nanocontacts>i& They have also been used successfully 
to explain the experimental current-voltage character- 
istics of molecular nanowires connecting non-magnetic 
metal electrodesiSi 3 ^ 3 ^ In this paper, we model the Fe 
nanocontacts as nanoclusters of Fe atoms connected to 
ideal leads that represent the source and drain electrodes. 
However we parameterize the electronic structure of the 
Fe clusters with the use of semi-empirical tight-binding 
parameters obtained from fitting the known band struc- 
ture of magnetic bulk Fe^ Never-the-less, as we explain 
in Section [n] our model also provides a satisfactory de- 
scription of the spin-resolved surface densities of states 
and of the local magnetic moments at the surfaces of Fe 
nanoclusters. Thus our model incorporates in an approx- 
imate way both the bulk and surface magnetic properties 
of Fe which together influence the SDT through a junc- 
tion of bulk Fe leads that come together at a nanocontact, 
as in experimental realizations of Fe atomic chains. 



This article is organized as follows: In Section ITT1 we 
describe our model of the geometry and electronic and 
magnetic structure of the Fe nanocontacts. In Sect ion ITU 
we summarize the formalism used in our spin-dependent 
transport calculations. In Section IIVI this formalism is 
applied to different arrangements of Fe atoms connecting 
the ferromagnetic nanoclusters: We start by considering 
a structure in which the positions of the atoms of the Fe 
nanocontacts and of the atomic bridge connecting them 
coincide with sites of an Fe crystal, however we consider 
a single atomic bridge, as distinct from the periodic array 
of nanocontacts treated in Ref. 0. The spin-dependent 
currents flowing through this junction are calculated uti- 
lizing initially a simple, voltage- independent model of the 
transmission probabilities of spin up and spin down elec- 
trons through the junction. A strong magneto-resistance 
is predicted. We then perform voltage-dependent calcu- 
lations of the transmission to examine the effects of finite 
bias on the spin-dependent predictions more closely and 
obtain a still strong, though relatively weaker, JMR. We 
then proceed to examine geometries in which the sepa- 
ration of nanoclusters is allowed to vary, stretching the 
atomic chain. To this end, we supplement the above 
tight-binding model of the electronic structure of the Fe 
clusters^ by introducing a position-dependent parame- 
terization for the electronic coupling between clusters, in 
the spirit of extended-Hiickel theory. We allow the sep- 
aration between the tips of the clusters to range from 
a bulk nearest-neighbor distance, in which case ballis- 
tic SDT is predicted, to a distance at which the atomic 
chain has been broken and vacuum-tunneling-like SDT 
occurs. The cross-over regime is examined and interest- 
ing phenomena are predicted, including spin-dependent 
transport resonances mediated by dangling bonds on the 
tip atoms. In the ballistic regime, unusual device char- 
acteristics such as negative JMR and spin polarization 
of the current are predicted. We then proceed to ap- 
ply our approach to the case of tunneling between ferro- 
magnetic nanoscale tips separated by vacuum and exam- 
ine the spin-dependent current and magneto-resistance. 
Physical mechanisms are presented that explain the pre- 
dicted effects. The I — V characteristics and magneto- 
resistance predictions presented in this paper should be 
experimentally accessible. Our conclusions are summa- 
rized in Section Ivl 



II. THE MODEL 

The system of interest consists of two bulk ferromag- 
netic metal electrodes that act as a source and drain for 
electrons, joined by a nanoscale junction of atoms of the 
same magnetic metal. We model this structure as a con- 
nected pair of nanoscale contacts of (100) body-centered 
cubic iron, each contact consisting of a 55 atom cluster 
built from 5x5, 4x4, 3x3, 2x2 layers of atoms, terminated 
with a single tip atom. Thus the atomic chain bridging 
the nanocontacts is a dimer consisting of the two tip Fe 
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atoms. 

The electronic structure of the Fe clusters is described 
by a tight-binding Hamiltonian using a non-orthogonal 
basis of s, p and d atomic orbitals, a total of 9 orbitals per 
atom for each spin orientation. The values of the Hamil- 
tonian matrix elements H^j and overlaps Sij = 
between orbitals i and j are, where possible, taken from 
Ref 36. These tight-binding parameters are based on fits 
to ab initio band structures calculated for Fe crystals 36 
and have previously been employed successfully to study 
magnetic multilayer systems^ and magnetic*^ and non- 
magnetic^ atomic chain nanoscale systems. However 
this parameter set is modified and/or extended appropri- 
ately (as discussed below and in Sections llV Al and llV B II 
respectively) when a magnetic domain wall is present or 
a finite bias voltage is applied across the junction or the 
positions of the atoms of the junction do not all coincide 
with sites of a perfect bcc Fe crystal lattice. 

The source and drain electrodes arc modeled as sets 
of one-dimensional ideal leads coupled to each cluster. 
These ideal leads are implemented as semi-infinite tight- 
binding chains with one orbital per site, one such lead 
being coupled to each (s,p and d, spin up and spin down) 
orbital of each atom of the two layers of each Fe cluster 
that are furthest from the junction. The site-orbital en- 
ergy of each ideal lead is matched to that of the Fe tight 
binding orbital to which the lead couples and the lead 
band widths are chosen large enough that all of the eigen- 
states of the leads are propagating modes in the energy 
range of interest. As well as mimicking macroscopic elec- 
trodes by supplying an ample electron flux to the system, 
this large number of ideal source and drain leads (369 per 
cluster for each spin orientation) has an effect similar to 
phase-randomizing Biittiker probes 38 in minimizing the 
influence of dimensional resonances due to the finite sizes 
of the Fe clusters employed in the model. Thus our re- 
sults do not change qualitatively if the parameters of the 
ideal leads are varied (within reasonable bounds) or an 
additional layer of Fe atoms is included in the clusters. 

We assume that each cluster (including the tip atom) 
and the semi-infinite leads connected to it form a single 
magnetic domain. Thus, if the magnetizations of the two 
nanocontacts are anti-parallel, an atomically thin, hard 
domain wall is present at the midpoint of the Fe dimer 
connecting the two clusters, in accordance with the cal- 
culation of Ref. |3^ and with previous work modeling 
SDT in Ni atomic contacts^ The tight-binding parame- 
ters developed in Ref. |3f| do not include the off-diagonal 
Hamiltonian or overlap matrix elements between atomic 
orbitals located on opposite sides of a magnetic domain 
wall. However the published off-diagonal matrix elements 
connecting spin up orbitals are not very different from the 
corresponding matrix elements connecting spin down or- 
bitals in the same domain.- 36 This is physically reasonable 
given the that spatial part of a spin up orbital should be 
similar to that of the corresponding spin down orbital. 
Thus in the present work, at zero bias, when geometri- 
cally applicable, we approximate the off-diagonal matrix 
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FIG. 1: (Color online) Calculated bulk (a-d) and surface (e) 
densities of states vs. energy (eV) at zero bias for spin up 
(solid curve u) and spin down (dotted curve d). The Fermi 
energy is located at OeV. (a) The densities of states due to 
d-electrons. (b) The densities of states due to p-electrons. 
(d) The densities of states due to s-electrons. (d) The total 
densities of states. The d orbital contribution dominates the 
total density of states except at the higher and lowest energies 
shown, (e) The surface contribution to the densities of states. 



elements connecting orbitals that are on opposite sides of 
an atomically thin domain wall and have the same spin 
orientation by the average of the spin up and spin down 
matrix elements connecting the corresponding orbitals in 
a single domain. As in in Ref. |3|| we also take all Hamil- 
tonian and overlap matrix elements connecting orbitals 
having opposite spin orientations to be zero, i.e., spin 
flips during electron transmission through the junction 
are not considered. 

To better understand the significance of the results of 



our SDT calculations it will be useful to compare them to 
particular features of the spin-resolved surface and bulk 
densities of states of the Fe electrodes and to further 
resolve the density of states into its s, p and d com- 
ponents. Since the atomic orbital basis that we use is 
non-orthogonal we carry out this resolution as follows: 
Consider the normalized eigenstate = X^ c mN) °f 
the Hamiltonian, where \i) is an atomic orbital. Express 
the norm of the eigenstate in terms of the contributions 
of all of the atomic orbitals of the system as 
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where Sij — To resolve the contribution of 

eigenstate \^k) to the density of states into its s, p 
and d orbital components we assign to each atomic or- 
bital i the Mulliken weight c* k i cu,i + Y^jjti( c k,i c k,jSi,j + 
c k j c k,iSj,i) /2. The partial and total bulk densities of 
states obtained in this way are shown in Fig. ^ ( a ~d) ■ 

In order to assess the applicability of our semi- 
empirical tight binding model (that is based on fits to 
bulk band structures^) to nanostructures that include 
surfaces, we calculated the local magnetic moment per 
atom in the surface atomic layer and in the interior of 
each of several large Fe clusters. A similar Mulliken anal- 
ysis to that described above was used to resolve the elec- 
tron probability distributions for individual eigenstates of 
the tight-binding Hamiltonian used in the present work 
into surface and interior contributions. We found the 
lower coordination number and lack of symmetry at the 
surface to result in enhanced magnetic moments at the 
surface, approximately 2.5/is per atom at the surface 
versus 2.2/zb in the bulk for our model, a result qualita- 
tively similar to that of ab initio surface calculations for 
Fe (2.25 in the bulk and 2.98 at the surface)i£ We note 
that Mulliken analyses arbitrarily assign half of the prob- 
ability contribution that is due to overlaps between two 
atoms in equation @ to each of the two atoms involved 
and therefore tend to underestimate differences between 
the electronic populations of adjacent atoms, often signif- 
icantly. Thus we regard the level of agreement between 
the present model and the ab initio surface calculations^ 
as satisfactory. The contributions from the surface layer 
of a representative Fe cluster to spin up and spin down 
densities of states obtained in this way are shown in Fig- 
ure ^e) and compare well with those reported in Ref. 

As further indication that our semi-empirical model is 
appropriate for describing also the nanocontact region 
(the two tip atoms) in addition to the rest of the con- 
tacts approximating bulk electrodes, we have employed 
the model to predict magnetic moments of actual Fe nan- 
oclusters of increasing size. The results of this calculation 
are compared to those of ab initio calculations on iden- 
tical nanoclusters 41 in Figure |3fa-d). It is evident that 
the predicted moments of our semi-empirical model com- 
pare well to those of the ab initio study for both the sur- 
face and interior atoms of the clusters, particularly when 
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FIG. 2: (Color online) Calculated magnetic moments per 
atom as a function of atom distance relative to the central 
atom of Fe nanoclusters. Solid lines correspond to the present 
study, dashed lines to the first principles calculations of Sipr 
et aim— All atoms are at bulk lattice positions and cluster size 
is increased as additional coordination shells are added, (a) 
15 atom Fe cluster. (b) 27 atom cluster, (c) 51 atom cluster, 
(d) 59 atom cluster. 



the cluster size resembles that of our nanoclusters used 
for the transport calculation (51 and 59 atom clusters 
compared to the 55 atom clusters we employ in our fol- 
lowing study). We consider the good agreement for the 
outermost atoms of the clusters as indication that our 
model adequately reproduces the magnetic properties of 
even atoms of low coordination number, such as the tip 
atoms of the clusters in our transport calculations. Ad- 
ditionally, as the nanocluster size is increased (to the 51 
and 59 atom sized clusters) , the magnetic moments of all 
of the interior atoms quickly approach that of the bulk. 
Therefore these results also provide further evidence that 
the Fe clusters employed in our transport calculations 
are large enough to model macroscopic magnetic Fe elec- 
trodes (with single atom tips) if our semi-empirical tight 
binding model is used to treat the electronic and mag- 
netic structures of these systems. 
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III. THEORY OF SPIN-DEPENDENT 
TRANSPORT 

When a bias V is applied between the ferromagnetic 
Fe electrodes an electric current I flows through the 
nanoscale junction described in Section [H] Landauer 
theoryS relates this current to the multi-channel proba- 
bility T for an electron to scatter from the source elec- 
trode to the drain via the junction, according to 

HV) = i J dET(E,V)\f(E,ns)-f(E,n D )] (4) 

where E is the energy of the electron, F(E, /i) is the 
equilibrium Fermi distribution and /j-s.d = Ep ± eV/2 
are the electro-chemical potentials of the source (S) and 
drain (D) electrodes in terms of the common Fermi en- 
ergy, EpSk The transmission probability in equation (JIJ 
is 

T(E,V)= ]T |^||t/Mw| a (5) 

where tp tS '- taiS is the transmission amplitude from a state 
of ideal lead a of the source electrode with spin s to a 
state of ideal lead (3 of the drain electrode with spin s' and 
v aiS and is/3. s > are the corresponding electron velocities. 
If spin flips during transmission of electrons through the 
junction are neglected as in the present work and the 
summations over spin in equation (JjjJ for T are restricted 
to a particular spin orientation then equation Q yields 
the current for that spin orientation, i.e., I up or Idown- 

We calculate T numerically by solving the Lippmann- 
Schwinger equation^ 

\V a ) = \<S>«) + G (E)W\y a ) (6) 

where Go (E) is the Green's function for the system with 
the coupling W between the ideal leads and Fe clusters 
switched off, |$q ) is the eigenstate of ideal lead a when it 
is decoupled from the Fe clusters, and \^ a ) is the corre- 
sponding eigenstate of the complete system with coupling 
W switched on. The transmission amplitudes tpy- i(x>s 
that appear in equation JSJ are obtained from the solu- 
tion of the Lippmann-Schwinger equation ©. 

IV. RESULTS 

A. Fe Nanocontacts Bridged by Fe Atoms in a bec 
Nearest-Neighbor Geometry 

We begin by considering the nanocontact geometry 
shown in the inset of Figure OUa). Here the two 55 atom 
clusters are placed so that the atom terminating the tip 
of each cluster is at a bulk nearest neighbor position rel- 
ative to the tip atom of the other cluster. Thus all atoms 
occupy positions that match sites of a body-centered cu- 
bic Fe bulk lattice so that a tight-binding Hamiltonian 
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FIG. 3: (Color online) (a) Transmission probabilities as a 
function of energy (eV) at zero bias for the contact geometry 
shown in the inset where all of the atoms are located at sites of 
a bulk Fe bec lattice. Spin up —* up (uu), down — > down (dd) 
transmission probabilities are for parallel magnetization of the 
two contacts; up — * down (ud) is for anti-parallel magnetiza- 
tion. The Fermi energy is located at OeV. (b) The current 
as a function of voltage for parallel (p) and anti-parallel (a) 
magnetizations calculated by integrating the zero bias trans- 
mission probabilities in (a). The JMR (j) is also shown, (c) 
The currents and JMR as calculated with an applied bias us- 
ing the linear voltage drop model, (d) The currents and JMR 
as calculated with an applied bias using the abrupt voltage 
drop model. 



and overlap matrix for the entire structure can be con- 
structed using the parameter set of Ref . [3(| Figure (3{a) 
shows the calculated transmission probabilities for the 
different spin configurations, at zero bias, for this geom- 
etry. The spin up — > up (uu) and down — > down (dd) 
transmission probabilities are for parallel magnetization 
of the two contacts; up — > down (ud) is for anti-parallel 
magnetization. The Fermi energy is at OeV. The trans- 
mission characteristics depend somewhat on the specific 
details of the chosen geometry, such as cluster size and 
shape. However, quite similar results are obtained when 
a 6x6 atom layer is added to each of the clusters. Cur- 
rent and magneto-resistance results are also robust to 
geometrical changes, such as addition of a 6x6 layer, or 
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deletion of the 5x5 layer that terminates the nanocluster, 
indicating that the clusters and ideal leads in our model 
adequately represent real (macroscopic) Fe leads. 

For this geometry, the calculated spin up — > up trans- 
mission, perhaps surprisingly for a transition metal, ex- 
hibits a fairly flat plateau close to unity near and above 
the Fermi energy. As for atomic chains constructed 
from some non-magnetic metals^ 7 -, this transmission near 
unity in the spin up case is not due to a single open 
conducting channel, but to the superposition of partially 
conducting channels. By systematically switching off the 
tip atom-to-tip atom coupling for different atomic orbital 
symmetries, we deduce that the up — > up transmission 
near the Fermi energy is due primarily to s and p elec- 
trons, even though as is seen in Figure^the spin up den- 
sity of states due to d electrons is comparable to those 
due to s and p electrons above the Fermi energy. Be- 
low the Fermi energy, the bulk and surface densities of 
states for up electrons are larger (primarily d-like), and 
the predicted transmission is also higher. 

In contrast to the spin up case, the spin down — > down 
transmission near the Fermi energy, not only shows a 
strong contribution from p-type electrons but a very 
strong d-electron component; s-electrons do not seem im- 
portant in this case. This also can be contrasted to the 
spin down densities of states very near the Fermi energy 
where that due to p-electrons is small and comparable to 
that for s-electrons. Yet the presence of many common 
features that can be observed in the calculated trans- 
mission and the surface and bulk densities of states does 
support the conclusion that the spin transport properties 
of this magnetic system are influenced by the densities of 
states. However, the transmission characteristics can not 
be determined from solely a knowledge of the densities of 
states; characteristics such as the near quantum conduc- 
tance above the Fermi energy in the spin up channel and 
the magnitudes of the transmission features can not be 
deduced from the densities of states alone. Those char- 
acteristics can be determined only from a full quantum 
transport calculation such as that presented here. 

Interestingly, as shown in Figure Ola), the calculated 
spin down — > down transmission is significantly larger 
than the spin up — > up transmission near the Fermi en- 
ergy and thus according to equations and I0J our 
theory predicts a negative spin polarization of the cur- 
rent through this junction at low bias for parallel mag- 
netization of the two contacts. Negative spin polar- 
ization of the current is a common feature of theories 
that attempt to estimate spin-dependent transport ef- 
fects in Ni and Co heterostructures simply from convo- 
lutions of the spin-dependent densities of states of the 
magnetic electrodes that are dominated by the contri- 
butions of d electrons. However, such approximations 
have been shown incorrect in magnetic tunnel junctions, 
as d-electron wave functions decay very quickly in the 
insulating b arr i e r^i2iLSL2i4iii2ii. This decay is further 
enhanced by the conservation of wavevector parallel to 
the interfac e) 45 ' 47 , and therefore the d-electrons that are 



responsible for the strongest features in the densities 
of states do not play an active role in transmission in 
such insulator mediated systems. In atomic chain sys- 
tems however, there is no lateral periodicity, and parallel 
wavevector conservation is broken. Also, the magnetic 
electrodes that we consider are in this case in physi- 
cal contact and so d-electrons do play an active role in 
transport. Never-the-less, a convolution of the bulk 
densities of states cannot account for the negative spin 
polarization of the current that we predict in the present 
system; it arises from a combination of effects due to the 
surface and bulk densities of states and our full quantum 
mechanical treatment of electron transmission through 
the junction. 

The transmission probability T(E, V) that enters the 
Landauer expression Q for the electric current / flow- 
ing through the nanocontact depends on the applied bias 
voltage V as well as the electron energy. The bias depen- 
dence of T depends on the potential profile through the 
nanocontact which is difficult to calculate from first prin- 
ciples since this is a non-equilibrium many-body prop- 
erty. However, appropriate heuristic models for the pro- 
file can yield accurate results for the current 4£ We adopt 
this approach here by comparing the results obtained for 
a variety of simple models. 

Figure IHIb) shows the current calculated by approx- 
imating T{E 1 V) with T(E, 0), an approximation com- 
monly used to study transport at low values of the bias. 
Here, and in all of the figures and discussion that follow, 
the current refers to the flow of electrons, or electron 
flux, as electrons are transported from the source to the 
drain electrode. In Figure EHb) the calculated current 
is higher for parallel magnetization than for anti-parallel 
magnetization and steadily increases with applied bias, 
characteristic of ballistic transmission. The calculated 
JMR, as defined by Eq. ^ is positive and larger at low 
bias, decaying to a lower value of about 0.15 at high bias. 

Figure |3Jc) and (d) show our results obtained us- 
ing explicitly bias-dependent transmission probabilities 
T(E, V) that were calculated for two different models of 
the potential profile of the nanocontact. In each case 
we assume that the entire voltage drop occurs across the 
narrowest constriction in the system, i.e, over the dimer. 
Thus all atoms (and ideal leads) to the source side of 
the dimer are assumed to be at a potential (f> = 
while atoms and leads to the drain side are at + 
For the results shown in Figure [2c) the potential is as- 
sumed to vary linearly through the region occupied by 
the dimen 4 ?i 5 9i 51 i 52 i 5 ? In Figure^d), for comparison, the 
potential is assumed to change abruptly from <p = ±y 
in the atomic layers adjacent to the dimer to = on 
the two atoms that constitute the dimer itself, a profile 
analogous to that proposed initially in theoretical work 
on ballistic semiconductor nanostructures 5 ^ and more re- 
cently adopted in modeling certain molecular wires^ 1 * 5 ^ 

In each case the electrostatic potential modifies the di- 
agonal matrix elements of the tight-binding Hamiltonian 
of the system which become Hi i = — e<j>i where if? ■ 
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is the tight-binding Hamiltonian matrix at zero bias and 
4>i is the electrostatic potential (j) at the site occupied by 
atomic orbital i. Because the tight binding basis that we 
use is non-orthogonal, the applied electrostatic poten- 
tial also modifies the non-diagonal Hamiltonian matrix 
elements^ Here we include this effect approximately in 
the form 

H iJ = H9.-eS iJ (4> i + 4,^/2. (7) 

The applied electrostatic potential breaks the symme- 
try between the left and right clusters in the inset in 
Figure EI a )- Symmetry breaking often results in weaker 
transmission probabilities in quantum transport, and this 
has qualitative implications for the present system: The 
symmetry breaking as energy levels are shifted apart 
manifests itself in a somewhat lower transmission and 
current for parallel magnetization of the contacts (es- 
pecially at higher bias) in Figure Et c ) and (d) than 
in Figure Et a ) where the effect of the applied bias on 
the transmission is neglected. However the net current 
for anti-parallel magnetization is much less sensitive to 
bias-related symmetry effects^ Because it selectively de- 
presses the current for parallel magnetization of the con- 
tacts, the bias-induced symmetry breaking results in a 
crossover with increasing bias from positive to negative 
values of the JMR in Figure Et c ) and (d) , an effect not 
found in the less realistic model of Figure Eta) where the 
effect of the bias on the electron transmission probability 
is neglected. 

The bias voltage at which we predict negative JMR to 
appear (~2V in Figure Et c ) an d ~3.5F in Figure Eld)) 
depends on the details of the potential profile across the 
junction where the two electrodes touch: In the linear 
voltage drop model the bias applied across the junction 
simultaneously shifts the energies of similar atomic or- 
bitals on the two tip atoms, bringing them closer together 
or further apart, depending on their relative spin orien- 
tations and whether the magnetizations of the contacts 
are parallel or antiparallel. Thus orbitals of the tip atoms 
are brought closer to or further from resonance with each 
other. On the other hand, since in the abrupt voltage 
drop model the applied bias affects the energies of all of 
the atomic orbitals of the contacts except those of the tip 
atoms, the corresponding resonant effect of the applied 
bias in this model is weaker: energies of orbitals within a 
cluster are shifted with respect to unshifted energy levels 
on the tip atoms (the energy levels of the rest of the clus- 
ters are not directly coupled to any other atoms in the 
opposing cluster). Therefore, the onset of negative JMR 
requires a higher applied bias in the abrupt voltage drop 
model. 



B. Spin-Dependent Transport for More General 
Junction Geometries 

It has been established experimentally that the con- 
ductance characteristics of atomic chains can be altered 



by stretching themS. We now investigate the corre- 
sponding dependence of the JMR in our magnetic sys- 
tem. 

1. Generalizing the Tight-Binding Model 

For the more general junction geometries, our Hamilto- 
nian matrix elements between atomic orbitals within each 
of the two Fe clusters will again be based, as discussed 
above, on the Fe ferromagnetic tight-binding parameters 
derived from ah initio bulk band structure calculations 36 . 
However a different approach is needed to obtain the 
non-diagonal Hamiltonian matrix elements that describe 
tip-to-tip coupling since the atoms involved no longer 
all fall on sites of a single bec Fe lattice. We estimate 
these matrix elements using an appropriate modifica- 
tion of extended Hiickel theory. Extended Huckel— is a 
semi-empirical tight-binding model from quantum chem- 
istry that provides a simple description of the electronic 
structures of a wide variety of molecules. It uses a non- 
orthogonal atomic orbital basis. The diagonal Hamilto- 
nian matrix elements are identified with the exper- 
imental ionization energies of the corresponding atomic 
orbitals i of isolated atoms. In the Wolfsberg-Helmholtz 
form of the modelSi the non-diagonal matrix elements 
are assumed to be 

K 

where Sij are the overlap matrix elements connecting 
orbitals i and j that are calculated by approximating 
the atomic orbitals i and j with linear combinations of 
Slater-type orbitals. K is a phenomenological parameter, 
usually chosen to be 1.75, in order to match experimental 
data. We estimate the non-diagonal Hamiltonian matrix 
elements that describe tip-to-tip coupling using equation 
(jHJ) and numerical values of the diagonal Hamiltonian ma- 
trix elements based on those given in Ref. y(| However, 
the values of the diagonal matrix elements given in Ref. 
l36l are not close to the ionization energies of the corre- 
sponding atomic orbitals but are defined up to an ar- 
bitrary additive constant. It is necessary to choose the 
value of this constant with care in order to obtain re- 
alistic results from equation iJSJ. We choose its value 
so that when the two tips are positioned in the geome- 
try studied in Section IIV Al , our model (including the 
tip-to-tip coupling matrix elements obtained using equa- 
tion JSj) reproduces the up — > up transmission near the 
Fermi energy obtained using only the bulk parameters 
from Ref. |3y. The value of the shift that we apply to 
the bulk electronic parameters is c = — lZ.hieVmL The 
same tip-to-tip coupling parameters H m ,n are then used 
for all spin configurations^ This is consistent with the 
fact that bulk up-up coupling parameters are very similar 
to down-down parameters^ 6 -. 

Thus in our calculations the overlaps Sij between 
Slater-type orbitals provide the distance and orientation- 
dependence of the tip-to-tip Hamiltonian matrix elements 
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required for the more general junction geometries (as in 
extended Hiickel theory), however, the Hamiltonian ma- 
trix elements are normalized so as to yield transport re- 
sults consistent with those obtained from a Hamiltonian 
matrix derived from ab initio calculations. 



2. Structural Considerations: Bulk and Relaxed Geometries 

The first geometry we consider using this hybrid model 
is shown in the inset of Figure 0Ja). The two (100) 
5x5, 4x4, 3x3, 2x2, and single tip atom clusters are now 
aligned with a common axis. The two tip atoms are again 
separated by a bulk nearest neighbor distance, 2.482A. 
We refer to this as the "linear geometry," and to the 
structure in the inset of Figure^a) as the "skewed geom- 
etry." Our geometry relaxations^ on simple model sys- 
tems involving the pair of 2x2 planes and the two atoms 
forming the dimer, indicate that the two tip atoms in 
this geometry are close to their stable positions. Sim- 
ilar relaxations that we carried out showed that when 
the clusters are pulled apart, some stretching occurs be- 
tween the tip atoms and their respective 2x2 layers, but 
most of the stretching occurs between the two tip atoms. 
When the tips are far apart representing vacuum tun- 
neling, our geometry relaxations show that, with respect 
to the rest of the cluster, the tip atom sits very near its 
bulk position. Therefore, we assume that as the junction 
is stretched, all of the stretching occurs between the two 
tip atoms, and so the atoms of each tip separately are 
located on their bulk Fe lattice positions. Therefore, in 
what follows, we again make use the bulk Fe ferromag- 
netic tight-binding parameters to describe the electronic 
structure within the clusters for all matrix elements in 
and use our modified Hiickel approach to describe 
the inter-cluster coupling^ 



3. SDT in Fe Nanocontacts Bridged by Fe Atoms in 
Contact in a Linear Geometry 

Our results for spin transport in the linear geometry 
with the tip atoms separated by a bulk nearest neigh- 
bor distance are shown in Figure 0] The spin up — > up 
transmission shown in Figurc^a) shows a similar plateau 
above Ep to that in the skewed geometry (Figure EJa)). 
This similarity is due to the dominance of the non- 
directional s-orbitals in both cases. Below the Fermi 
energy, the transmissions exhibit significant differences, 
due to the importance of the very directional d-orbitals in 
that energy range. For the two different geometries, spin 
down — > down and up — * down transmissions are similar 
in a broad sense, such as overall magnitude, but display 
many differences on a finer scale due to the importance of 
d orbitals in those transmissions. Figure^Jb) shows that 
the current through this geometry, as calculated from the 
voltage-independent transmissions, is smaller than that 
through the skewed geometry using bulk parameters to 
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FIG. 4: (Color online) (a) The transmission probabilities as a 
function of energy (eV) at zero bias for the spin up — > up (uu), 
down — > down (dd), and up — > down (ud) configurations for 
the Fe atoms at the nearest neighbor separation, arranged 
in the linear geometry (inset). The Fermi energy is located 
at OeV. (b) The current as a function of voltage for parallel 
(p) and anti-parallel (a) magnetizations in the zero bias ap- 
proximation for the transmission probabilities. The JMR (j) 
is also shown, (c) The currents and JMR as calculated with 
an applied bias using the linear voltage drop model, (d) The 
currents and JMR as calculated with an applied bias using 
the abrupt voltage drop model. 



describe the coupling. However, with second nearest- 
neighbor couplings across the junction turned off in the 
skewed geometry, the current in the parallel magnetiza- 
tion configuration is in very close agreement for the two 
geometries suggesting that the lower current can be at- 
tributed to the lack of second nearest neighbor coupling 
between the two contacts in the linear geometry. In the 
anti-parallel magnetization configuration, the current is 
substantially lower in the linear geometry than in the 
skewed geometry. This reflects the lower spin up — > down 
transmission in the the linear geometry which results 
partly from weaker d coupling: Many of the overlaps in- 
volving d orbitals are zero in the linear geometry. The 
substantially lower anti-parallel current results in a sig- 
nificantly stronger JMR. 
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As shown in Figure^c) the current for parallel magne- 
tization and the magneto-resistance, as calculated from 
the voltage-dependent transmission in the linear volt- 
age drop model, are again significantly lower than for 
the voltage-independent transmission model for the lin- 
ear geometry (at higher bias) because the Hamiltonian is 
no longer symmetric. The current is slightly weaker for 
the linear geometry than for the skewed geometry and a 
larger JMR is predicted than for the skewed geometry. A 
negative JMR is again predicted, but this time it mani- 
fests near a bias of 3V. As was discussed for the skewed 
geometry, in the linear voltage drop model the applied 
bias simultaneously moves the atomic orbitals of the tip 
atoms closer or further apart in energy. Thus the linear 
voltage drop model provides a relatively strong mecha- 
nism for those energy levels to move closer to or further 
from resonance with each other. However, in the present 
case of the linear geometry, there is no second neighbor 
coupling between a tip atom and atoms of the opposite 
contact and so the resonant effect is weaker than in the 
previous case of the skewed geometry. The potential pro- 
file with abrupt voltage drops is also less conducive to 
negative JMR here (as in the skewed geometry) thus a 
negative JMR does not occur within the voltage range 
in Figure Eld). It should also be noted that the predic- 
tions regarding I — V characteristics and JMR in the two 
geometries are qualitatively similar for the linear volt- 
age drop model which, according to density functional 
calculations,*- 9 may be a more accurate approximation for 
all-metal systems than the abrupt voltage drop model. 
However, the predicted results are qualitatively similar 
for both models of the potential profile in both systems, 
demonstrating the robustness of our method; it is rea- 
sonable to expect the results of a fully self-consistent cal- 
culation of the potential profile to fall between these two 
model profiles. Again, all of these results are also reason- 
ably robust to the addition of a 6x6 atom layer to each 
of the clusters. 



4- Dependence of the Spin Transport on the Separation 
Between Tip Atoms and Dangling Bond Formation 
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FIG. 5: (Color online) Spin dependent transport in the lin- 
ear tip geometry vs. separation of the tip atoms, (a) The 
JMR as a function of tip separation calculated using the lin- 
ear voltage drop model and applied biases of 0.5V and 1.0V. 
(b) The spin polarization of the current as a function of tip 
separation for parallel magnetization of the contacts, (c) The 
spin polarization of the current as a function of tip separa- 
tion for anti-parallel contact magnetizations, (d) The energy 
eigenvalue associated with the spin down dangling bond state 
relative to the chemical potential of the drain electrode (OeV) 
as a function of tip separation for parallel magnetization of 
the contacts. 



Figure EJVc) shows the calculated dependence of the 
magnetoresistance and spin polarization of the current 
(as defined by Eqs. Q and (J2J, respectively) on the 
separation between the two tip atoms as the junction 
is stretched for two values (0.5 and 1.0V) of the bias 
voltage, for the linear voltage drop model. 

As the separation between tip atoms is increased, the 
JMR (Figure Ufa)) initially decreases rapidly then in- 
creases to a local maximum at tip separations near 3.5 
and 4.A then resumes its decrease turning weakly negative 
and, for the 0.5V bias, increases again at large separa- 
tions. The spin polarization of the current for parallel 
magnetization of the contacts (Figure Etb)) is negative 
at small separations, i.e., the spin down (minority spin) 
current predominates, and rises initially with increasing 



separation almost to zero near 3A. It then passes through 
a minimum near 3.5 and 4 A (where the JMR shows a 
maximum) before resuming its rise and becoming posi- 
tive at large separations. The local maxima in the JMR 
and minima in the spin polarization of the current near 
3.5 and 4 A are due to a pronounced resonance in the dom- 
inant down — > down transmission that appears near the 
Fermi energy for parallel magnetization of the contacts 
at such separations. This transmission resonance persists 
when a 6x6 Fe layer is added to each nanocontact, thus 
it does not appear to be a dimensional resonance due to 
the finite size of our Fe clusters. 

By examining the contributions of the individual elec- 
tronic eigenstates of the Hamiltonian of the coupled Fe 
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clusters to the down —> down transmission (through the 
Green's function Go in the Lippmann-Schwinger equa- 
tion ©) we have identified the particular eigenstate that 
is responsible for this transmission resonance. Our Mul- 
liken analysis revealed that for tip separations corre- 
sponding to the appearance of the down — > down trans- 
mission resonance, a significant portion of this eigenstate 
resides on the tip atom of the drain contact and has 
d(x 2 — y 2 ) orbital symmetry there; the d(x 2 — y 2 ) or- 
bital has a lobe oriented along the s-axis towards the tip 
atom on the other cluster. Since the d(x 2 — y 2 ) orbital is 
involved in tip-to-tip bonding and the amplitude of the 
eigenstate of interest on the tip atom increases as the 
tip separation is increased, we attribute the appearance 
of this eigenstate state to bond-breaking between the tip 
atoms and the associated formation of a dangling bond. 

The energy E of the dangling bond eigenstate relative 
to the electrochemical potential of the drain is plotted 
in Figure |5{d) as a function of tip separation for the two 
values of applied bias considered in Fig. GJa-c). The dan- 
gling bond state is the dominant feature in the Mulliken 
spectra of the tip atom within a broad window about the 
Fermi energy and this criterion is used to identify the 
state for the different values of the tip separation and ap- 
plied bias. At zero temperature, only states within the 
limits of integration corresponding to the applied bias 
contribute to the predicted current (see Eq. These 
limits are determined by the electrochemical potentials 
of the source (0.5eV and l.OeV in Figure^d) for the two 
values of applied bias) and the drain (OeV). As shown in 
the figure, as the tip separation is increased, the energy 
of the dangling bond state shifts into the window of inte- 
gration, above the electrochemical potential of the drain 
and therefore begins to contribute to the current. (Note 
that this happens at separations close to those where the 
JMR in Fig[SJa) first begins to rise signaling the onset of 
the transport resonance). Simultaneously, the Mulliken 
weight of the state located on the drain tip increases as 
the separation is increased to a value of about 10 percent, 
a sizeable portion of the total probability distribution of 
the eigenstate given that a total of 110 atoms make up 
the Fe clusters in our model. Therefore, as the separa- 
tion is increased from its smallest distance of 2.482A, we 
attribute the initial decrease in JMR to the stretching 
of the bond between the two tip atoms. The following 
increase to the local maximum in JMR is attributed to 
the formation of a dangling bond as the tips are fur- 
ther pulled apart, leading to the resonant feature in the 
down — > down transmission. Once the dangling bond has 
been formed, its energy and Mulliken weight are roughly 
constant, and the JMR resumes its decrease with increas- 
ing separation. 

As was discussed in Section IIV Al s, p and d elec- 
trons all play a significant role in transport in this sys- 
tem for small tip separations. However their contribu- 
tions decay differently as the tip separation is increased 
and the tunneling regime is entered. Since the valence d- 
electrons have a lower site energy than the other electrons 



(Ed < E s ,E p )r^ they decay more rapidly resulting in a 
less significant contribution from d-electrons to transport 
as the tip separation is increased. 

For an applied bias of 0.5V and near a tip separa- 
tion of about 4.5A, the contribution to tunneling from 
d-electrons is roughly equal to that of s-electrons. Here 
Tup^up, mostly due to s-electrons, and T down ^ down , 
mostly due to d-electrons, are roughly equal and the spin 
polarization of the current in the parallel magnetization 
case (Fig. [S] (b)) is roughly zero, defining a cross-over 
regime. This can be viewed as the cross-over from bal- 
listic transmission to a tunneling- like transport. At the 
cross-over, transmission from spin up (mostly s and p) 
to spin down (mostly d) for antiparallel magnetization of 
the contacts is roughly the same as that for the parallel 
spin configurations, and the resulting JMR (Fig. El (a), 
0.5V) is very small, or even slightly negative. Since this 
small negative magneto-resistance appears at a relatively 
small applied bias, it should be accessible experimentally 
in these systems. 

At larger separations, beyond the cross-over regime, 
the current is dominated by the transport of s and p elec- 
trons and T up ^ up makes the largest contribution. There- 
fore the spin polarization of the current in the parallel 
magnetization case (Fig. (b)) is positive, and grow- 
ing. Since transmission involving down (mostly d) elec- 
trons is rapidly decreasing, T down ->down and T up 

— >down 

are small (due to band mismatch the anti-parallel trans- 
mission in general is smaller than even the decreasing 
down — > down transmission) and, as can be seen in Fig[5] 

(a) in the 0.5V case, the JMR also begins to slowly grow. 
At higher biases (for example at 1.0 V bias in Fig [SI a) and 

(b) ) more energy levels are sampled, increasing the im- 
portance of the d states and the cross-over doesn't occur 
until larger separations. Therefore, this model predicts 
that when d states are important, the current will be 
dominated by minority electrons, and when d states are 
negligible, current due to majority electrons will domi- 
nate. This is quite analogous to the predictions made 
by Mathon on periodic systems involving Co and tun- 
neling gaps of varying widths, where the fast decrease in 
d-electron transport accounts for a rapid reversal in sign 
of the spin-polarization^. 

Since the application of bias breaks the symmetry of 
the Hamiltonian, T uv ^ down ^ T down ^ uv and a non-zero 
spin polarization of the current is predicted in the anti- 
parallel magnetization case as well. This effect is entirely 
due the voltage drop across the junction and so is ex- 
pected to be unique to magnetic transport systems that 
can sustain a significant potential drop across the junc- 
tion, i.e, atomic nanocontacts and tunnel junctions. As 
shown in Figure [S^c), the spin polarization of the cur- 
rent for antiparallel magnetization of the contacts (as 
defined by Eq. [2Jl is negative over a wide range of tip 
separations due to the dominance of the up — > down 
transmission under applied bias as a result of shifting 
energy levels. Even at relatively low bias, there is more 
source and drain energy level overlap in this configura- 
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FIG. 6: (Color online) (a) The transmission probabilities as a 
function of energy (eV) at zero bias for the spin up — > up (uu), 
down — > down (dd), and up — > down (ud) configurations for 
a 7A tip separation representative of the vacuum tunneling 
regime. The Fermi energy is located at OeV. (b) The cur- 
rent as a function of voltage for parallel (p) and anti-parallel 
(a) magnetizations as calculated by integrating the zero bias 
transmission probabilities in (a). The JMR (j) is also shown, 
(c) The currents and JMR as calculated with an applied bias 
using the linear voltage drop model. 



tion than in the down —> up configuration where levels 
are shifted apart by the applied bias. At higher bias, this 
is even more the case, resulting in a much more negative 
spin polarization of the current at 1.01^. Never-the-less, 
the spin polarization does increase for large tip separa- 
tions as the d-state transport becomes less significant. 
The non-monotonic behavior of the spin polarization of 
the current with tip separation for separations below 4A 
in Fig[SJc) is due to resonant transport associated with 
dangling bond formation similar to that for the case of 
parallel magnetization of the contacts that is discussed 
in detail above. 



5. Vacuum Tunneling Between Fe Nanocontacts 

When the separation between tip atoms is very large 
(vacuum tunneling), the d-clcctron states have decayed 
across the gap between the contacts and do not signif- 
icantly contribute to the transmission. This situation 
is easier to analyze as the couplings between tip atoms 
involving s and p x (all the rest are negligible) behave 
simply. Figure EJa) shows the calculated transmission 
for the linear geometry with a tip atom separation of 



7A at zero bias. The currents calculated using the zero 
bias approximation for T, shown in Figure [S^b), display 
step-like features due to peaks in the transmission. The 
zero-bias approximation for T also predicts a strong JMR 
effect, which shows broad peaks where the transmission 
resonances are encountered. 

Figure 0c) shows the I — V characteristics calculated 
using the linear voltage drop model for T(E, V). The 
application of bias and shifting of energy levels strongly 
damps the local maxima in the zero bias transmissions 
and the step-like behavior seen in Figure E{b) is now 
less pronounced. The current increases with bias, resem- 
bling the accepted behavior of tunneling through vac- 
uum. At the application of about 2V, the up — > down 
transmission becomes the strongest scattering channel, so 
the anti-parallel current becomes stronger than the par- 
allel current. This results in a relatively strong, negative 
JMR, although the down — > down transmission regains 
its dominance and positive JMR returns at higher bias. 
Thus the shifting of different transmission maxima again 
results in non-monotonic behavior of the JMR and neg- 
ative JMR for some values of the applied bias. However, 
there is no sustained JMR reversal such as is predicted 
when d-electron states contribute to transport. 



V. CONCLUSIONS 

We have presented a microscopic quantum theory of 
spin-dependent transport across iron nanoscale junctions 
bridged by chains of Fe atoms, based on the Landauer ap- 
proach to transport, semi-empirical tight-binding Hamil- 
tonians and Lippmann-Schwinger and Green's function 
scattering techniques. We first applied bulk ferromag- 
netic tight-binding parameters (that were shown to also 
provide a satisfactory description of surface properties) 
to study ballistic transport between a pair of Fe con- 
tacts connected by two Fe atoms in a nearest neighbor 
geometry. We presented theoretical predictions for the 
current-voltage characteristics of this system for paral- 
lel and anti-parallel magnetizations of the contacts and 
predicted that negative spin polarization of the current 
should occur at low bias. We also predicted that the 
junction magnetoresistance of this system should switch 
sign from positive to negative with increasing bias. Next, 
we extended our model so as to allow us to study spin- 
dependent transport for more general tip geometries and 
presented predictions of the junction magnetoresistance 
and spin polarization of the electric current through Fe 
nanoscale junctions as a function of separation between 
the tip atoms of the nanocontacts from nearest neigh- 
bor distances to the vacuum tunneling regime. Char- 
acteristic trends emerging in those transport predictions 
as the separation is varied were associated with decay 
rates of different orbitals and with the breaking of bonds 
and associated dangling bond formation resulting in spin- 
dependent transmission resonances. We also presented a 
systematic physical interpretation of our predictions that 
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in many cases are qualitatively different from those that 
rely on density of states considerations alone, and also 
different from the behavior of SDT in thin film junctions 
where lattice periodicity in directions orthogonal to that 
of the current flow plays a key role. Many of our pre- 
dictions apply in the regime of low bias that should be 
accessible with present day experimental techniques. 
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